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The distribution function of the local density of states is 
computed exactly for the Wigner-Dyson ensemble of random 
Hamiltonians. In the absence of time-reversal symmetry, pre- 
cise agreement is obtained with the "supersymmetry" theory 
by Efetov and Prigodin of the NMR lineshape in disordered 
metal particles. Upon breaking time-reversal symmetry, the 
variance of the Knight shift in the smallest particles is reduced 
by a universal factor of 2/3. 

PACS numbers: 76.60. Cq, 71.20.Ad, 73.20.Dx 



A new quantum size effect in small metal particles has 
been predicted by Efetov and Prigodin ||l|]. They com- 
puted the spectral hneshape for nuclear magnetic res- 
onance (NMR) and found that the resonance becomes 
very broad upon decreasing the temperature and parti- 
cle size, due to large fluctuations in the Knight shift at 
different points in the sample. Similar results for a dis- 
ordered wire have been reported by Mirlin and Fyodorov 
1^, who extended earlier work on one-dimensional con- 
ductors |^-||] . Recent NMR spectroscopy on a monodis- 
perse array of nanometer-size Pt particles is in striking 
agreement with the theory |Q. The essential difference 
with quantum size effects on thermodynamic properties 
1^ is that NMR in a metal measures the local density 
of states p{E,t) = Y.n^{^ - ^ri)\^ n{r)\'^ , and thus de- 
pends both on the energy levels £"„ and the wavefunc- 
tions 5'n(r) of the valence electrons. The sensitivity to 
the local density of states occurs because the frequency of 
microwave absorption by a nucleus is shifted as a result 
of the Fermi contact interaction between the nuclear spin 
and the electron spin. (This is known as the Knight shift 
l^.) The frequency shift 6uji for a nucleus at position 
is linearly proportional to p(Ep,ri) (with Ep the Fermi 
energy). Different nuclei have different frequency shifts, 
which is observed as a broadening of the resonance in an 
ensemble of particles. To determine the broadening one 
has to consider the combined effect of particle-to-particle 
fluctuations in the energy spectrum and spatial fluctua- 
tions of the wavefunctions within the disordered particles. 
The strong spectral correlations make the theory consid- 
erably more complex than for an ideal spherical particle 
00. 

Efetov and Prigodin computed the fluctuations in the 
local density of states from a microscopic model of a dis- 
ordered metal particle with elastic impurity scattering 
in a strong magnetic field, by mapping it onto a super- 
symmetric non-linear sigma-model. The traditional ap- 



proach 1^ to quantum size effects in metal particles is 
based on random-matrix theory. Following Gor'kov and 
Eliashberg , it is assumed that the Hamiltonian 7i of 
an irregularly shaped metal particle is a random N x N 
Hermitian matrix, with the Gaussian distribution 



PiH) = Cexp(-cTrH2) 



(1) 



originally introduced by Wigner and Dyson for the spec- 
trum of a heavy nucleus |l3| . The coefficient c determines 
the mean level spacing A (which in the limit iV — > oo is 
uniform in the bulk of the spectrum) , and C is a normal- 
ization constant. In the past, random-matrix theory has 
been applied to quantum size effects on thermodynamic 
properties of metal particles ||^,[T^ , in agreement with mi- 
croscopic theories p5|-p^ . These applications involve the 
distribution of the set of eigenvalues {En} of H, which 
follows from P{'H) on integrating out the eigenvectors 
{^„}. In contrast, the NMR lineshape depends on the 
joint distribution of the E's and ^f's. Problems of this 
type have not previously been tackled by random-matrix 
theory. 

It is the purpose of this paper to show how the anoma- 
lous broadening of the NMR line shape can be obtained 
directly from the Wigner-Dyson distribution (|l|) , without 
any further assumption. In the absence of time-reversal 
symmetry we recover precisely the results of Ref . . Ex- 
periments on nanometer-size particles are typically per- 
formed in the presence of time-reversal symmetry. (The 
authors of Ref. estimate that to break time-reversal 
symmetry in their system would require magnetic fields 
of the order of 1000 T, two orders of magnitude greater 
than the experimental fields.) Random-matrix theory is 
particularly suited to investigate the dependence of the 
fluctuations on fundamental symmetries of the Hamilto- 
nian. A celebrated example is the reduction by a factor 
of 1/2 of the variance of the universal conductance fluc- 
tuations, upon breaking time-reversal symmetry [p7|-p^ . 
We will show that the variance of the Knight shift has a 
different reduction factor of 2/3, provided the level spac- 
ing is much greater than both the level broadening and 
the temperature. 

Let us first reformulate the problem of the NMR line- 
shape in the framework of random-matrix theory. The 
intensity /(w) of the resonance at frequency lo is given by 
the distribution P{p) of the local density of states upon 
rescaling. 



aP[p = b{uj - ujq)), 



(2) 
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with microscopic parameters a,b,LUQ. The distribution 
P{p) is defined by 



(3) 



where the average (• ■ •) is a spatial average over the total 
volume occupied by the particles. The local density of 
states p{E, r) is given by 



p{E,r)^j:j{E-En)\'fn{r) 
f(E)^{^l2^){E^ + W)-\ 



(4) 
(5) 



where 7 is the broadening of the levels due to tunneling 
into the medium in which the particles are imbedded. 
Eqs. (|^) and (^) assume that 7 is greater than the tem- 
perature, the electronic Zeeman energy, and the spin- 
orbit scattering rate. For completely isolated particles, 
even-odd electron number effects play a role ||^, which 
are not considered here. (These effects are expected to be 
relatively unimportant in the metal-cluster compounds of 
current experimental interest [Q.) 

For an iV-dimensional Hamiltonian Ti, the continuous 
variable r is replaced by the index m = 1, 2, ... A^, and 
\E'„(r)p becomes {N/V)\Umn\^ i with U the unitary ma- 
trix that diagonalizes 7i and V the volume of a parti- 
cle. In the absence of time-reversal symmetry, U varies 
over the full unitary group. This is relevant to NMR 
for very strong magnetic fields and not too small parti- 
cles, and is the case considered in Ref. [Q. If the flux 
penetrating a particle is much less [ pO| than h/e, then 
time-reversal symmetry is not broken and U is restricted 
to the orthogonal group. The orthogonal and unitary en- 
sembles are characterized by the index (3 — 1,2, which 
counts how many real numbers Umn,q il — 1, ■ • ■ /^) define 
the matrix element Umn- There exists a third symmetry 
class, characterized by /3 = 4 and U a symplectic matrix, 
which describes systems with time-reversal symmetry in 
the presence of strong spin-orbit scattering All 
our /3-dependent formulas for the local density of states 
apply also to the symplectic ensemble, however, the ap- 
plication to NMR requires a modification of the theory 
because spin and charge density are no longer directly 
related. 

For each of the random-matrix ensembles, the average 
in Eq. (^) can be written as an integration over eigenval- 
ues and eigenvectors 

P{p) = JdEi---J dEN J due exp ^ E^j x 



i<j 



Y[\E,-E,f5ip-^Y.f(^-^\U-' 



V 



(6) 



Here, and in what follows, we choose Ep as the zero 
of energy. The Jacobian Yl^^j \Ei — Ej\^ introduces /3- 
dependent correlations between the eigenvalues, in the 
form of level repulsion. The eigenvectors are uncorre- 
lated with the eigenvalues, and distributed uniformly 
with measure dU. 



The variance of the Knight shift requires the first two 
moments of P{p)- The general formula for the p-th mo- 
ment is 



(7) 



w = (N/vr ({j:j{e^)\u^„\^Y 

The first moment evaluates trivially to 

p={VAr'=Po. (8) 
To evaluate the second moment we use the formula 



U„ 



\Umn'\ ) 



N{l3N + 2)' 

In the limit N — + 00, at constant V and Z\, we find 



(9) 



pVPo - 1 + (1 + 2//5)Z\ / dEf{E)^ 

'J —00 

/OO pOO 
dE dE'f{E)f{E')T2[E-E'). (10) 
-00 J —00 

The two-level cluster function 

UE - E') = (E.^/(S - E^E' Ej)) (11) 

is known |2^. For j3 — 2 one has 

T2{E)^{'kE)-'^s\t?{ttE/A). (12) 

The expressions for /3 = 1 , 4 are a little more complicated 
1^ . The asymptotic behavior of Eq. ([l0| ) for Z\ ^ 7 is 
obtained from MmE^Q A^T2{E) ~ 1, hence 

1*00 

7lpl = {l + 2lf5)A dEf{E) + 0{j/A). (13) 



In the opposite regime Z\ ^ 7 one may approximate 
AT2{E) ~S{E), hence 

OA 

p'/pI - 1 + / dEf{E) + 0{Ahf. (14) 



From Eq. (|10|) one readily computes the variance 
YaiK/k = ~(?IpI - 1 of the Knight shift. The result 
is plotted in Fig. 1 for /3 = 1, 2 (solid curves). The small 
and large-Z\ asymptotes (dashed) are both linear, but 
with different slopes: 



,- 2A f ^ if Zi<7, 



(15) 



We have checked that the values for /? = 2 agree with 
Ref. [|]. The transition /3=l^/3 = 2on increasing 
the magnetic field is signaled by a reduction of Var KjK 
by a factor of 2/3 for Z\ > 7 and 1/2 for Z\ < 7. The 
reduction by 1/2 is the same as for the variance of the 
universal conductance fluctuations (UCF) in mesoscopic 
metals. In that case the broadening of the levels is always 
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much greater than the level spacing. (Their ratio is the 
conductance in units of e^/ft., which is 3> 1 in a metal.) 
The reduction by 2/3 has no analogue for UCF. 

So far we have assumed that the level broadening 7 is 
much greater than the thermal energy k^T. At higher 
temperatures, the fmiction f{E) in Eq. (^o|) is to be re- 
placed by the convolution /t(-B) of the Lorcntzian ^ 
with the derivative of the Fermi function. Its Fourier 
transform frik) = / dE e''''^ fT{E) is 

/T(fc) = 7rfcBTfc[exp(i7|fc|)sinh(7rfcBrfc)]"\ (16) 

The variance of the Knight shift for 7 ^ k-oT is plotted 
also in Fig. 1 (dotted curves). The asymptotic formulas 
for small and large A become 



A 



kp,T V knT 



7 



if Z\ < 7 + k^T, 



^ ifZi»7 + fcBr, 



$(s) 



Jo 



(17) 



We conclude that the reduction factor associated with 
f3 = 1 P = 2 is universal, independent of the relative 
magnitude of temperature and level broadening. This 
is relevant for experiments, which are typically in the 
regime that k^T and 7 are of comparable magnitude [Q. 

We now turn to the complete distribution P{p), given 
by Eq. (^ . It is convenient to work with the (dimension- 
less) Laplace transform 



F{s) 



dp exp{-sp/po)P{p), 



(18) 



and recover P{p) at the end by inverting the transform. 
First, we average over the eigenvectors. It is known 
p3| , p3[ that, to leading order in 1/A^, the PN compo- 
nents u„in,q {n — 1,2, ... N; q — 1, . . . /?) of a single row 
of U are independently distributed Gaussian variables 
with zero mean and variance 1 / PN. Carrying out the 
Gaussian integrations, we find 



(19) 
(20) 



g{E,s) = l-[l + {2sA/P)f{E)]-P/'. 



The remaining average over the eigenvalues can be car- 
ried out using the method of orthogonal polynomials pl| . 
This method works for any /3, but is simplest for the case 
/3 = 2. In that case the function g is a Lorentzian in E, 

g{E,s) = {sjA/2n){E' + \r'r\ (21) 

with r = 7(1 + 2sZ\/7r7)i/2. The largc-iV limit of Eq. 
( p^ ) for P — 2 is given by the Fredholm determinant 

n^) = nr=i(i-A«(s)), (22) 

where An is an eigenvalue of the integral equation 



/oo 
dE' g{E', s)tI'\e - E')4>{E') = \4>{E). 
-00 



Fourier transformation gives 



2tt 



-1x1 A 



dk' g{k — k' , s)(j){k') — A0(fc), 



5(A:,s) = (,s7Zi/r)exp(-ir|fc|). 



(23) 

(24) 
(25) 



To evaluate the Fredholm determinant of Eq. (24), we 
discretize k £ (-7r/Z\, vr/Zi) as kn = {n/A){-l + 2n/M), 
n = 1, 2, . . . M, and then take the limit M ^ 00 p41: 



F(s) = lim Det 



S7 



: exp 



TrFln 



MT V MA 
(cosha + i(7/r + r/7)sinha). (26) 



Inversion of the Laplace transform yields finally 

P{p) ^ i(a/2^)i/2a;-3/2 exp[-ia(a; + x'^)] 

X (cosha -I- i(x -I- — a~^) sinha). (27) 

Here x = p/po and a = nj/A. Equation (p^) is precisely 
the distribution of Efetov and Prigodin |^J7^ 

This solves completely the problem for /3 = 2 and zero 
temperature. For /3 = 1,4 and T 7^ the distribution 
function can still be written as a Fredholm determinant, 
which then has to be evaluated numerically. The Wigner- 
Dyson distribution (|l|) can only describe the pure sym- 
metry classes {P — 1,2, or 4). The transition between 
symmetry classes might be studied by means of an ex- 
tension known as Dyson's Brownian motion model pl| . 
We leave these problems for future work. 

In summary, we have derived the result of Efetov 
and Prigodin Q for the NMR lineshape in the absence 
of time-reversal symmetry from the single assumption 
that the Hamiltonian of the particle is a member of the 
Wigner-Dyson ensemble of random-matrix theory. A 2/3 
reduction of the variance of the Knight shift in the small- 
est particles has been predicted to occur upon breaking 
time-reversal symmetry. 

I am indebted to K. B. Efetov for suggesting this prob- 
lem as a challenge for random- matrix theory. B. Rejaei 
helped me to compute the determinant in Eq. (26). Dis- 
cussions with H. B. Brom on the experimental aspects 
have been most helpful. This work was supported by the 
Dutch Science Foundation NWO/FOM. 
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FIG. 1. Dependence of the variance of the Knight shift K 
(normalized by its average K) on the level spacing A, with 
and without time-reversal symmetry (/3 = 1 and 2) . The solid 
curves are computed from Eq. (|lo|), for the case that the level 
broadening 7 is much greater than the thermal energy ksT. 
The dotted curves are for the opposite case 7 <^ k^T. The 
dashed lines are the asymptotes ( p^ ) and ( [l4| ) for small and 
large A. (The factor 6/7r in the definition of T* is chosen such 
that the solid and dotted curves have the same asymptotes.) 
Breaking of time-reversal symmetry reduces Var if by a factor 
of 1/2 and 2/3 for small and large A, respectively. These 
factors are universal, independent of the relative magnitude 
of 7 and fceT. 
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